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It has long been thought that superconductivity breaks down even at zero temperature in lower- 
dimensional systems due to enhanced quantum phase fluctuations. In ID wires, these fluctuations 
have been described in terms of "quantum phase-slip" (QPS): tunneling of the superconducting order 
parameter between states whose relative phase differs by ±2n. Many deviations from conventional 
bulk superconducting behavior have been observed in ultra-narrow superconducting nanowires over 
the last several decades which have been identified with QPS, and a number of theories have been 
developed which appear to be supported by at least some of the observations. However, other 
observations in many cases point to contradictory conclusions or cannot be explained by existing 
theories, such that a unified understanding of the nature of quantum phase slip and its relationship 
to the various observations has yet to be achieved. In this paper we present a new model for QPS 
which builds on an idea recently postulated by Mooij and Nazarov [Nature Physics 2, 169 (2006)]: 
that flux-charge duality, a classical symmetry of Maxwell's equations, can be used to relate QPS 
to the well-known effect of Josephson tunneling. We hypothesize that this duality holds at the 
microscopic level, and based on this we build a theory which predicts the fundamental QPS energy 
scale identified by Mooij and Nazarov and observed in recent experiments, and which can explain a 
variety of the other observations made on these quasi- ID systems. The results of our theory, and the 
corresponding interpretations of observations which they generate, are qualitatively different from 
previous treatments. They may also provide a conceptual link to quantum phase fluctuations in 2D 
superconductors. 

PACS numbers: Valid PACS appear here 
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Phase fluctuations are one of the important ways that 
mean-field theories of superconductivity break down in 
lower-dimensional systems, such as lattice planes in high- 
Tc superconductors [l], , thin films > & n d narrow 
wires [6j. In quasi- ID wires, they come in the form of 
"phase slips" : topological excitations in which the phase 
difference <j> between the wire's ends changes by ±2ir. 
Such a process was first suggested by Little [7| , in which 
thermal energy causes a short segment of the wire to fluc- 
tuate into the normal state, allowing <f> to "slip" by one 
cycle. This was later treated by Langer, Ambegaokar, 
McCumber, and Halperin (LAMH) [8[ and others [9J as 
a source of finite resistance just below the critical temper- 
ature Tc, and subsequent experiments [l(| on ^0.2-0.5 
/im-diameter crystalline Sn "whiskers" validated these 
ideas. 

In 1986, by analogy to macroscopic quantum tunneling 



in Josephson junctions (JJs) [ll|, Mooij and co-workers 
proposed that for narrow enough wires, a quantum phase 
slip (QPS) process might exist, in which the superconduc- 
tor tunnels between states whose phase differs by ±27r 
12]. Shortly thereafter, using lithographically defined, 
~ 50 nm-wide In wires, Giordano measured finite resis- 
tance that persisted much farther below Tc than LAMH 
predicted [1J], exhibiting a crossover from LAMH scal- 
ing with T — Tc to a slower temperature dependence. 
This was interpreted as direct observation of a crossover 
from thermal to quantum phase fluctuations. A pioneer- 
ing theory was later developed by Golubev, Zaikin, and 
co-workers (GZ) [l4| which had some success in support- 
ing this intuition with microscopic physics, and many 



researchers have since built on these ideas [27], [28|, |47J , 
and used them to try to connect newer observations to 

qps |, EMI. 

However, in other recent experiments on Pb [l8l . [2l| 
and MoGe [f| [23| nanowires < 10 nm wide, the 
anomalous low-T resistance crossover was in many cases 
not observed at all. This is difficult to reconcile with 
theories that associate the crossover directly with QPS, 
since it should depend exponentially on the wire's cross- 
sectional area, and it has since been suggested that the 
observed crossover may actually be associated only with 
granularity [IH[l|l[24[ and/or inhomogeneity [25[ of the 
wires rather than QPS. Even more striking was an appar- 
ent complete destruction of superconductivity as T — > 
in some of these nanowires with a normal-state resistance 
R„ >Rq = h/4e 2 0, H3, [23jj2^1. Although theories exist 
which predict insulating |27H29| or metallic [Tl l30l. l3lj 
states in ID as T —> 0, it remains unclear whether any 
can explain a T = critical point at R n ~ Rq. Thus, 
a unified understanding of the various observations and 
their connection to QPS has not yet been achieved. 

In 2006, Mooij and Nazarov (MN) [32[ made what may 
turn out to be a conceptual leap forward: they postu- 
lated that a classical symmetry known as flux-charge du- 
ality [33M36] ] can be applied to connect QPS with Joseph- 
son tunneling, the well-known process in which Cooper 
pairs penetrate through a thin insulating barrier sep- 
arating two superconducting electrodes, and establish 
phase coherence between them. Based on this idea, MN 
posited the existence of a phase slip potential energy 
Ups(q) = —Escosq, dual to the Josephson potential 
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Uj(4>) — —Ejcos(j). Here, <f> and q are known as the 
phase and quasicharge, Ej is the well-known Josephson 
energy, and E$ is a new energy scale for QPS, which MN 
left as an input parameter. From this elegant hypothe- 
sis, they generated a phenomenology of QPS dual to that 
of JJs, including a dual set of classical nonlinear equa- 
tions for q [371 ] . and a dual class of circuits involving ID 
superconducting nanowires, what they called "phase-slip 
junctions" (PSJs) [3^, [H, [&| . Very recently, in fact, two 
of these circuits have been demonstrated experimentally 
EH j including what appear to be direct observations 
of the phase slip energy Es, providing the most direct 
evidence yet seen for QPS in continuous wires fl9j . 

In this work, we describe a new theory which builds 
on and extends MN's postulate, and results in a unified 
model that appears to answer many of the outstanding 
questions. We take as a hypothesis what we call micro- 
scopic flux-charge duality between QPS and Josephson 
tunneling, an extension of MN's postulate which allows 
it to be applied to the calculation of the phase slip en- 
ergy Es itself. This is contained in section U below. Our 
result for Es is qualitatively different from previous the- 
ories, in two important ways. First, it centrally involves 
the dielectric permittivity e inside the superconductor, a 
quantity which does not appear in this way in any previ- 
ous theory for QPS. In our model the permittivity plays 
the role of an effective mass for "fluxons", fictitious dy- 
namical particles dual to Cooper pairs, whose motion 
"through" a ID wire corresponds to quantum phase slip, 
just as Cooper pair motion through an insulating barrier 
corresponds to Josephson tunneling. Second, our result 
implies that QPS involves fluctuations of the phase only, 
and not the amplitude, of the superconducting order pa- 
rameter. By contrast, previous theories such as that of 
GZ assume QPS involves not only phase fluctuations but 
also necessarily fluctuation of the order parameter am- 
plitude to zero; this inevitably results in a more compli- 
cated, dissipative process. Although we do not attempt 
to prove the validity of the duality hypothesis in this pa- 
per, the predictions it generates for the phase slip energy 
agr ee q uantitatively with recent direct measurements of 
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In section UH we build on these results to construct a 
distributed, nonlinear transmission line model of a super- 
conducting nanowire, going beyond the lumped-element 
model of MN. We show that in the presence of QPS, its 
dynamical equations for quasi-classical phase evolution in 
one spatial and one time dimension (1+1D) can be cast 
into a form identical to the Maxwell-London equations 
in two spatial dimensions (2D), and from this establish 
a direct analogy between electric flux penetration into 
a superconductor in 1+1D and magnetic flux penetra- 
tion in 2D. In the process, we introduce a new length 
scale called the electric penetration length Xe, which is 
dual to the Josephson penetration depth in a long J J. We 
then use this analogy to predict quasi-classical phase ex- 
citations in 1+1D which we call type I and type II phase 
slips, which are direct electric analogs of a mixed state 



in a type I superconductor and the magnetic vortex in a 
type II superconductor, respectively. These type I and II 
phase slips are secondary macroscopic quantum processes 
|34| , in the sense that they arise as a collective effect out 
of the "primary" , microscopic QPS process, just as Bloch 
oscillations arise out of Josephson tunneling [34], H(| • 

In section IIII1 we introduce a simple model for the 
nanowire's electromagnetic environment, similar to that 
used previously for JJs [54|. We use this in conjunction 
with our transmission line model to provide a new in- 
terpretation of the "quantum temperatures" observed by 
Bezryadin [Hj], which were taken as direct evidence for 
QPS. We provide for the first time, to our knowledge, 
a quantitative explanation of the quantum temperature 
data, in terms of zero-point quantum fluctuations in the 
LRC oscillator formed by the nanowire and its environ- 
ment. Interestingly, in our model these observations de- 
pend only indirectly on the presence of QPS and very 
little on its actual magnitude Es- An important element 
of our explanation is the effect of a low environmental 
impedance at high frequencies, which provides damping 
for quantum phase fluctuations, and makes a description 
in terms of a quasi-classical phase appropriate. Related 
ideas were discussed previously by MN (331, and also in 
the context of JJs 

Closer to Tc, we show that our model also provides 
an alternative explanation for the crossover observed by 
Giordano [13j. Rather than a transition from thermal 
to quantum fluctuations as the temperature is reduced, 
it can be described in our model purely by a tempera- 
ture dependence of the energy barrier for thermal phase 
fluctuations, which itself undergoes a crossover when 
becomes longer than the wire length I. This new inter- 
pretation also provides an explanation for the mysteri- 
ous fact that Giordano-type crossovers are not observed 
in many wires with much smaller cross-sectional areas 
A and coherence lengths £ (for example, Bezryadin's 
MoGe wires have ^50 times smaller A and ~10 times 
shorter £ than Giordano's indium wires), even though 
QPS should increase exponentially as either of these pa- 
rameters are made smaller. In our model, the absence of 
these crossovers for some very narrow wires occurs natu- 
rally if stays shorter than I over the relevant temper- 
ature range. As the temperature approaches Tc, and the 
energy gap goes to zero, our model continuously connects 
to LAMH theory as one of its limiting cases. It also pro- 
vides an intuitive explanation for the dynamics that drive 
the wire's phase over the so-called "saddle-point" energy 
barrier [tHS, HH in terms of charging and discharging of 
the wire's kinetic capacitance. 

Finally, in section IIV1 we construct a possible expla- 
nation for the observed destruction of superconductivity 
when R n > Rq [23|], based on virtual type II phase slip- 
anti phase slip pairs as a fundamental quantum excita- 
tion. Unlike previous attempts to explain the apparent 
insulating behavior in terms of a dissipative phase tran- 
sition [Til . [2?i [28j | . we hypothesize instead a disorder- 
driven superconductor- insulator transition (SIT), closely 
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connected to phenomena observed in some quasi-2D su- 
perconductors 3 5, 42]. 
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FLUX-CHARGE DUALITY AND QUANTUM 
PHASE SLIP 



Flux-charge duality is a symmetry of Maxwell's equa- 
tions which can be derived purely classically. For 
lumped-element circuits, it manifests itself in the invari- 
ance of the equations of motion under the duality trans- 
formation of fig. [TJa). In continuous systems, it can be 
expressed in terms of the "quasicharge" Qs (defined for a 
surface S) and "fluxoid" $r (defined for a curve T) [43| : 



<3e = / atia-da), j q = j + — (i) 

E = -VK-f (2, 



$r = j> dt(E • ds) 



Figures [TJb) and (c) expand on this duality, and show 
how equations[T]and[2]can both be interpreted as a sum of 
contributions from "free" and "bound" quasicharge and 
fluxoid current densities Jq and J$: 
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where pq and B^- are the free charge and flux densities, 
and vq and v$ their velocities. Using the London gauge 
A = — AJq for a superconductor (where the kinetic in- 
ductivity is A = p$\ 2 with A the magnetic penetration 
depth) and D = eE for an insulator, yields: 



superconductor: A- 



dt 



insulator: 
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1 dt 



(5) 



(6) 



In a superconductor, free charge moves ballistically ac- 
cording to eq. [SJ which is London's first equation; in an 
insulator, free fluxoid can be viewed as moving ballisti- 
cally according to eq. [51 Maxwell's equation for the dis- 
placement current. Therefore, at the classical level of the 
Maxwell and London equations, superconductors and in- 
sulators are dual to each other. 

We now arrive at the proposed duality between a JJ 
and a PSJ, first suggested by MN, though here we have 
arrived at it in a different way. As shown in fig. [51 a JJ 




FIG. 1: Flux-charge duality, (a) duality transformation for 
lumped-element circuits; (b) and (c) the continuous case, (b) 
The free current density Jq = pqvq is the motion of free 
charge density pq at a velocity vq, through a surface £. The 
bound current density Jq = dD/dt is the displacement cur- 
rent density on E. (c) An illustrative example of "free" flux 
density B/, using a permanent magnet moving at velocity v$ 
relative to the stationary curve F, with J^ = v$ x Bf . In this 
construction, Jj, • ds is precisely the flux per unit time pass- 
ing through a segment ds. The bound fluxoid current density 
J|, = —dA/dt is associated with time- varying currents flow- 
ing along F, and the associated voltages from Faraday's law. 
Although the case of a moving magnet is somewhat artificial, 
any electric field in a medium can be broken into these two 
components: J^, associated with fields generated by bound 
charges, and j|, associated with induced emfs from time vary- 
ing currents (free charges) flowing in the medium. 



consists of two superconducting islands of Cooper pairs 
separated by an insulating potential barrier, while a PSJ 
is two insulating "islands" of fluxoid quanta (henceforth 
referred to as "fluxons" ) separated by a superconducting 
potential barrier. If we place the surface £ inside the 
insulating barrier of a JJ [Fig. G2a)] with junction ca- 
pacitance Cj, and the curve T inside a superconducting 
nanowire [Fig. [3jb)] of kinetic inductance Lk (we neglect 
the geometric inductance), we have: 



Qjj 

&PSJ 



n2e + CjV 




(7) 
(8) 



Here, n is the number of Cooper pairs that have passed 
through the JJ barrier, and CjV is the charge on the 
capacitance Cj of the junction. Similarly, L^I is the 
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total fluxoid of the nanowire associated with a current 
/, and m is the number of fluxons that have "passed 
through" the wire ($o = h/2e is the superconducting 
flux quantum). For thick enough superconducting wires, 
the only way for m to be nonzero is if some part of the 
wire was in the normal state at some time, as occurs in an 
LAMH phase slip over a length of wire ~ £, the coherence 
length. These events are dissipative, produce a measur- 
able voltage pulse, and can be associated with passage of 
a fluxon through the null in the superconducting order 
parameter at a localized, measurable position and time. 
By contrast, we take QPS to be coherent, delocalized 
fluxon tunneling through the entire length of wire, such 
that the order parameter remains nonzero and uniform 
in magnitude, no dissipation occurs, and no information 
about where the phase-slip occurred exists. Such a pro- 
cess would be the exact dual of Josephson tunneling. It is 
also qualitatively different from previous theories fl4Ll27l ]. 
in which a QPS event consists of a "core" region, inside 
which the order parameter is suppressed to zero, and out- 
side of which the phase dynamics are those of a linear 
transmission line. In such theories, the topological com- 
ponent (i.e. the actual phase "slip" itself) is presumed 
to be completely contained inside the normal-state core. 
Note that our discussion above based on the Maxwell 
and London equations does not prove the correctness of 
our alternative picture for QPS. Such a proof would re- 
quire a full microscopic theory dual to that used for JJs 
51], which we do not provide here. In this work, we 
assume such a theory could be formulated, but confine 
ourselves to exploring the consequences of our hypothesis 
and comparing them to experimental observations. The 
favorable results of this comparison, discussed in detail 
below, provide a strong indication that our hypothesis 
may be correct. 

We now seek to calculate the phase-slip energy Es , the 
central parameter for QPS. This quantity has been iden- 
tified [H Hi] with the QPS "rate" estimated by Giordano 
fl3l | , and later calculated by GZ (3] ■ In more recent the- 
ories it appears in the so called "QPS fugacity" and is 
either taken as the GZ result [13, HH H3 or left as an 
unknown [H, H3] . It is also the essential input parame- 
ter to the work of MN, the dual of the Josephson energy 
Ej for a J J |32j. We use the hypothesis of microscopic 
flux-charge duality to calculate it, using the well-known 
physics of Josephson tunneling as a guide. 

The Josephson effect can be understood intuitively 
by using a "macroscopic wavefunction" IS [49j for the 
superconducting state inside the barrier [50j . shown as 
iPq(x) in fig. Etc). Because of the finite Cooper pair mass 
(associated with the kinetic inductivity A), it is energeti- 
cally unfavorable for i/)q(x) to suddenly drop to zero upon 
entering the insulator (acquiring fourier components with 
large momenta), since the resulting zero-point kinetic en- 
ergy would far exceed the potential energy saved by ex- 
cluding the Cooper pairs from the insulator. Instead, the 
total energy is minimized by allowing ipq (x) to penetrate 
(tunnel) into the insulator over a distance Kq , such that 




FIG. 2: Flux-charge duality and quantum phase slip. Su- 
perconductor is shown in blue, insulator in red, and tunnel 
barriers are enclosed by black dashed lines, (a) an isolated 
JJ consisting of an insulating tunnel barrier between two 
superconducting islands with phases cj>\ and <j>2- The rele- 
vant canonically conjugate quantum variables are the relative 
phase 4> — <t>\ ~ <i>2, and the number of Cooper pairs which 
have tunneled through the barrier n. (b) an "isolated" PSJ 
consisting of a superconducting nanowire tunnel barrier be- 
tween two loops (fluxoid islands) with loop quasicharges q\ 
and q2- Here the canonically conjugate quantum variables 
are the nanowire's quasicharge q — q\ — q^ and the number of 
fluxoid quanta that have tunneled through the wire m. Note 
that while <j>\ , cf>2 , qi , q-i for the isolated islands in (a) and (b) 
cannot be well-defined, quasi-classical quantities, their differ- 
ences <j> and q can be. In (c) and (d) we add an electro- 
magnetic environment, in terms of an admittance Y env or an 
impedance Z env , with which one of the islands can exchange 
quanta. This results effectively in a single island coupled via a 
tunnel barrier to a reservoir. When Y env — > (c) becomes an 
ideal charge qubit with well-defined n and when Z env — > (d) 
becomes an ideal phase-slip qubit [44|, |45|| with well-defined 

the average kinetic and potential energies are equal. 

If microscopic flux charge duality holds between QPS 
and Josephson tunneling, then QPS can be understood in 
a dual manner as penetration of fluxons into a supercon- 
ductor [fig. [3jd)]. As described below, fluxons are ficti- 
tious particles with nonzero effective mass, each of whose 
motion through the wire corresponds to a phase slip. To 
illustrate the nature of this effective mass, and the result- 
ing quantum zero-point fluctuations which make fluxon 
tunneling possible in our model, we can use the example 
of fig- Die). The electrostatic energy density of this sys- 
tem can be re-expressed as a kinetic energy density for 
fluxons moving at velocity 

e|E| 2 V 2 
i Lp<$> 

where V = eE x B is the electromagnetic momentum 
density, E = v$ x B/ (in this example J = 0) , and 
p$ = e|B| 2 is an effective mass density. Based on this, 
and the identification of moving flux with time- varying 
$ [fig. HJc)], we identify the zero-point fluctuations re- 
sponsible for QPS with an effective mass for the fluxons 
associated with the permittivity e (dual to the kinetic 
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FIG. 3: Microscopic flux-charge duality between Josephson 
tunneling and quantum phase slip. Superconductor is shown 
in blue, insulator in red. (a) a JJ with the surface E inside 
the insulating barrier. In this case Q = C.jV + n(2e) where 
n is the number of Cooper pairs that have tunneled through 
the barrier; (b) a PSJ with T inside the nanowire, such that 
$ = Lkl + m&o where m is the number of fluxoid quanta 
( "fluxons" ) that have tunneled through the wire. Note that 
fluxons are one-dimensional objects [4g], shown for illustra- 
tive purposes only as rings, (c) Illustration of the macroscopic 
wavefunction rpQ for Cooper pairs inside the JJ barrier of 
thickness d, which can be understood as a superposition of 
two amplitudes: one with phase <f> penetrating inward in the 
+x direction (shown in red), and one with zero phase pene- 
trating inward in the —x direction (blue), (d) The proposed 
macroscopic wavefunction ip$ for fluxons defined inside the 
superconducting wire of radius R, which is a superposition 
of an amplitude with quasicharge q = Q/Qo penetrating in- 
ward in the —r direction (shown in red), and an amplitude 
with quasicharge zero penetrating outward in the +r direction 
(blue). These amplitudes are written in terms of the modi- 
fied Bessel functions Ko,Io- The effective mass of the Cooper 
pairs is associated with the kinetic inductivity A, which must 
be "charged up" for them to move. The effective mass of the 
fluxons is associated with the real part of the permittivity 
e inside the superconductor, which must be charged up for 
fluxoid to "move through" (for the phase to fluctuate) . 



inductivity A) . If the fluxon density were to drop to zero 
suddenly just inside the superconductor, this would re- 
quire fourier components of a fluxon wavefunction with 
large momenta (electric fields) and a large zero-point en- 
ergy cost. Just as with Cooper pairs in a JJ, the total 
energy is reduced when ip$, (r) extends into the supercon- 
ductor over a length shown in fig. EJd) . Note that 
the fluxons whose zero point motion is described by (r) 
are fictitious; they are formally dual to Cooper pairs, and 
just like Cooper pairs they have an effective mass, but in 
the end they are only a conceptual tool used to describe 
the macroscopic, coherent quantum phase fluctuations 
inside the ID superconductor. 



To evaluate the decay coefficient k$, we take V'$( r ) to 
(approximately) satisfy the Schrodinger-like equation: 



2eA 



U C A 



i/>$(r)=0, r<R 



(10) 



where A = irR 2 is the cross-sectional area of the wire, 
Uc is the condensation energy per volume of the super- 
conductor, and we take Q = ~iQoRd r (with the reduced 
charge quantum Q = 2e/27r), which is dual to the flux- 
oid operator for a J J barrier of thickness d: $ = —i<&add x 
(with the reduced flux quantum $o = <fro/2ir). In writing 
eq. [10] we have neglected interactions between fluxons, 
just as interactions between Cooper pairs are neglected 
inside a J J barrier [48]. The electric field can now be 
obtained from the ID fluxon probability current density 
associated with the 2D wavefunction u(r) = y/rijj^(r), 
thus (El: 



E=4 = 



Qo Rlm[u*(r)d r u(r)} 
71 27r|u(i?)| 2 



= E c smq (11) 



where q = Q/Qo, and Eq is the critical electric field, dual 
to the critical current density Jc for a JJ. The equality on 
the right can be deduced from ip<&( r ) i n n g-ISId); however, 
it is assumed to hold independent of the wire's cross- 
sectional shape (based on duality [Hj]). The phase-slip 
energy Eg = QqEqI for a wire of length I can now be 
written: 



Es s» lQo\ exp 



(12) 



where the approximate equality holds when the exponent 
is large [52| (which is also likely implied by our having 
neglected interactions between fluxons). Equation [12] in- 
dicates that materials with low Uc and e are desirable for 
strong QPS |53j . We defer comparison of this result with 
experiments to section IIII1 however, until after we have 
defined the lumped-element limit in which it is directly 
applicable. 

Note that eq. [12] also provides an intuitive explanation 
for the kinetic capacitance Cko = Qo/EqI discussed by 
MN (32|, which is dual to Josephson inductance. To see 
this, we write it in terms of a permittivity eko = Ckol/A: 



ffco 



K$R 



K&R > 1 



(13) 



This suggests that kinetic capacitance arises directly 
from the usual permittivity e (more precisely, its real 
part) inside the wire, magnified by a larger effective 
fluxon mass associated with tunneling through the po- 
tential barrier (the factor in parentheses). Although we 
do not treat the Uc — > (k^R — > 0) limit here, we ex- 
pect that in this limit efe — ¥ e and the fluxons see a simple 
dielectric rod (in this limit, purely classical flux-charge 
duality applies exactly). 
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FIG. 4: Model of a PSJ and its electromagnetic environ- 
ment. We treat the distributed ID system as a transmis- 
sion line, where in addition to the usual series (kinetic) in- 
ductance CkAx and shunt capacitance C±Ax per length Ax, 
we add a QPS amplitude that can be viewed as a (nonlin- 
ear) series kinetic capacitance Ck/Ax — Cko/Axsinq, with 
Cko = Qo/Ec (in units of Farads x length) and potential en- 
ergy (Qo/Cko)Ax cos q ^3^- This model is the exact dual 
of a long Josephson junction, with Ck dual to the Joseph- 
son inductance (Henry x length) and the polarization charge 
on C± dual to magnetic flux through the geometric induc- 
tance of the J J barrier (b) Electromagnetic environment 
model, following [H^]. The PSJ is indicated with the nota- 
tion of ref. ,32] (shaded box) with a kinetic inductor Lk in 
series with a QPS element, described by (lumped-element) 
energies E L = §l/2L k and E s (dual to E c = (2ef/2Cj 
and Ej for a JJ). The shunt capacitance C s h is associated 
with on-chip connections to the sample, Rdc with an exter- 
nal current source, and Rhf with the high-frequency char- 
acteristic impedance looking out of the package in which 
the sample is mounted. In typical experiments, we have: 
Rdc 3> Rq ^>Rhf- We assume also that ui v RhfChf ^> 1 
such that Z env (uj p ) w Rhf <C -Rq. 



now included free charge Qf (the imaginary part of the 
permittivity) as well. 

We can now describe fluctuations of t he sy stem shown 
in fig. d] using the partition function [lil |27|: 



Z = J X>*cxp[-5(*)/ft] 



where in the continuum limit Ax 
(imaginary-time) action is: 



(14) 

the Euclidean 



S 



dr / dx 



c, 



cosq 



in 



(15) 



and p = 1/ksT, t = it. The functional integral in cq.[T4l 
is carried out over paths VP in x, t which begin (at r = 0) 
and end (at r = h/3) in the uniform superconducting 
state (about which we consider fluctuations) . We assume 
that Gaussian fluctuations can be factorized out in eq. [14] 
(such that they simply renormalize the bare parameter 
values in S |27|), leaving only topologically nontrivial 
paths to be evaluated. 

To identify these paths, we note that the Wick rotation 
t — > —it inherent in eq. [15] can be viewed as a transfor- 
mation of a 1+1D classical dynamics problem into a 2D 
statics problem, and define: 



y = c s rt 



(16) 



where c s = y/l/ 'CkC± is the Mooij-Schon velocity [56| for 
electromagnetic wave propagation on the wire. We then 
obtain the following 2D equation of "motion" 57]: 



II. DISTRIBUTED QUANTUM PHASE SLIP 
JUNCTIONS 

The phase slip energy Es [eq. [T^] can be associated 
with a "bare" quantum phase slip element [32|, in the 
same manner that the Josephson coupling energy Ej for 
a JJ includes only the bare Josephson coupling and not 
yet the effects of the junction barrier's specific capaci- 
tance and geometric inductance. The dual quantities for 
a distributed nanowire are its series kinetic inductance 
per length Ck and geometric shunt capacitance per length 
to ground C±, respectively. We incorporate these into a 
transmission line model, shown in fig-SJa), in which QPS 
for a length Ax is shown as a nonlinear series kinetic ca- 
pacitance Ck/Ax = Cfeo/Axsing, with Cko = Qo/Ec- 
Note that in calculating Es and Ck above, we included 
only the effect of bound charges Q h (the real part of the 
permittivity), which describes the dielectric energy cost 
of the electric fields generated by phase fluctuations. By 
adding a kinetic inductance CkAx in series with each 
"bare" phase slip element Ck/Ax in fig. 0|a), we have 



A|V^ + sing = (17) 

where V^j, = x<9 x + yd y , we have used I/c s = Qad y q 
and C±V = —Qod x q, and defined an electric penetration 
length: 




and QPS plasma frequency ui p = \J\] CkCko |32|| . which 
are dual to the Josephson penetration depth and plasma 
frequency. We can now rewrite eq. [T7] for the dynamics 
of the system in fig. EJa) thus: 

Z L V xy x d = j (19) 
E 

Vxy X j = (20) 

with Z L = y^Ck/C^ the linear impedance of the wire, 
and the definitions: 
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FIG. 5: Type II phase slip in a ID superconductor. In 1+1D, a normal core of size £ is surrounded by circulating "currents" 
j [eqs. 1211 123| plotted in (a). A possible curve a for the line integral of eq. [52]is shown in red. Panels (b)-(f) show the order 
parameter along the wire as a phasor, at different times. On the leading edge of the vortex, current begins to flow over a 
length ~ A_b in the +x direction (b); this begins to charge up Ck, producing a gradient in j y — ZlP±- Once the total fiuxoid 
across the wire comes close to $o/2 (c), the order parameter can evolve continuously to a state (d) in which the current is zero 
and there is a null at the center of the vortex of length ~ £. This is the core of the type II phase slip and is related to the 
saddle-point solution encountered in long weak links [64| and for LAMH phase slips 0, Hj]. At this point the order parameter 
"passes through" the x-axis, and the supercurrent reverses (e). The current in the — x direction then discharges the kinetic 
capacitance as it ramps down to zero (f). This "static" solution in our x, y coordinates corresponds to an instanton-like solution 
[Hl> E3i HH m x,t. We emphasize that this is a quasi-classical solution for fiuxoid evolution by ±$o, which is conceptually 
distinct from QPS in the same way that Bloch oscillation in a JJ [3(| is distinct from Josephson tunneling. 



E = Ez 

d = Qz = e fc AE 

j = V xy ^ = £ k Ix + Z L p ± y (21) 

Here, p±_ = C\V is the charge per length on Cj_, and 
continuity requires that V xy ■ j = (note that the z 
direction is fictitious here, and is only introduced to use 
the curl notation in eas. 11912(1)) . Eauations[T^land[2^1have 
an identical form to Ampere's law and London's second 
equation in 2D, which would govern the penetration of 
perpendicular magnetic field into a thin superconducting 
film [481 ] : however, these equations govern the penetration 
of longitudinal electric field into a superconducting wire 
in 1+1D. In fact, it turns out that nearly all of the well- 
known magnetic results for the 2D thin film case have 
analogs here, including an electric Meissner effect, and 
type I and type II electric flux penetration [52j . 

We focus on the (electric) type II limit where A_e ^> £, 
which in most cases is appropriate for T — > (where Uc 
is largest and £ is smallest). We seek an electric vortex 
solution to eqs. ITT)ll2"T1 in which a normal core of size ~ £ 
in x, y is surrounded by circulating screening currents j 
[cq. [21] extending out to a radius ~ Xe- In order to 
include only closed paths in eqs. Q31 and 1151 we impose 
the condition (analogous to fiuxoid quantization in the 
2D magnetic case): 

* j • ds + / — • da = ±$ (22) 

J cr J a Cs 



where a is a closed curve in the xy plane which contains 
the core and bounds the surface a [fig. Ufa)]. This con- 
dition means that the fiuxoid evolves by $o (~^o) over 
the course of the event, which we call a "type II phase 
slip" (anti-phase slip). Using eqs. IT9"ll2"2l and assuming 
£fc ~ e/co far from the core (analogous to A( J) « A(0) for 
a magnetic vortex [Hf), we obtain [fig. [5]: 

j(p) « ±p-K 1 (■£-)$, P»e (23) 



where K is the modified Bessel function, p = ^/ x 1 + y 2 , 
and we assume h(3 ^> X E / c s = uj^ 1 . We can also deduce 
an interaction between type II phase slips separated by 
5p = \pi - p 2 \: 

where the sign is negative for a phase slip-anti phase slip 
pair. Both of these results are analogous to their mag- 
netic vortex counterparts [48j |. and this can be exploited 
to understand their implications. 



III. COMPARISON WITH EXPERIMENTS 

Before connecting our work to experiments, we must 
first include the effect of the electromagnetic environment 
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FIG. 6: Type II phase-slip in a low-impedance environment, (a) and (b) show the two lowest-energy Type II phase slip "lattices" 
for a constant voltage V — $ across a short wire (I < A_e), when R env (u} p ) <C Rq, Zl (These are analogous to Abrikosov 
lattices of magnetic vortices in superconducting strips of finite width). Image phase slips are shown with dashed lines, and 
result in currents that are nearly uniform along the wire, except near the cores of the type II phase slips (blue circles), (c) two 
lowest PS J energy bands Uo(&) and C/i(<3?) (which are exactly dual to the quasicharge bands of a JJ in a high-Z environment 
[34r - [36l l ). Inductive parabolae with E — El(&/&o — m) 2 are degenerate at E — El/4, where an avoided crossing occurs due 
to Es [3^. Eil. |45|] . If Es — > 0, the wire is simply an inductance Lk with energy El(&/<&q) (dashed black line). The current 
distribution shown in (a) for a short wire corresponds to adiabatic evolution along the lowest band shown by the dashed arrow 
in (c), which is dual to Bloch oscillation of a J J [34j - r36| . If the traversal remains adiabatic, the dynamics are insensitive to 
the magnitude of Es- Note that AE is smaller than the LAMH barrier by a factor of ~ £/Z, due to work done by the source 
[52| . and for Es <SC El is essentially the energy to charge up Lk with $ = $o/2 fr38l| . (d) shows the current-phase relation 
for the nanowire, which is nearly the same as for a long superconducting weak link [64[, but with additional avoided crossings 
associated with QPS, which produce a switching current I aw < Ic into a voltage state, (e) lowest two calculated energy bands 
Uo(Ib,&) and Ui(Ib,&) for wire SI of ref. [^ at lb —2 fiA 52]. (f) expanded view of the residual potential well in Uo{Ib, <&)• 
Fluctuations of the Lk — Rh f — C a h circuit produced by the wire and its environment can cause the phase particle to escape 
from this well even when there is still a potential barrier, at which point a voltage appears [ll|, UK • (g) calculated quantum 
temperature Tq [521 ] for wires Sl-5 of ref. [22I ] (based on ref. J6(J) vs. the observed Tq at T = 0.37^ [32]. With the exception 
of wire S3, the agreement is excellent, where the free parameters were (i) -R e „„(aj p ), for which the linear fit gives 115f2, and 
(ii) Cah, to which the results are insensitive as long as the system is overdamped (RhfC s h < VLkC s h), which is true here for 
Csh ;$ 10 fF. (h) Solid symbols are calculated values [52| of law plotted vs. the observed switching currents [67| for wires Sl-5 
in ref. 22] . The predictions are derived from eq. 1251 assuming that switching occurs when the potential well depth is equal 
to the observed To [531 . Two adjustable parameters were used (i) a constant a used to determine Lk from R n according to: 
Lk = ahRn/kBTc [4J]; and (ii) e. With a = 0.21 and e = 25eo, we find good agreement except again for wire S3. An estimate 
of a = 0.14 appears in ref. [4J| and e is not known. 



|58| . Our model is shown in fig. Hfb) , whose most impor- 
tant feature is a low, resistive impedance Rh / J$ y Ho/eo 
at the phase slip plasma frequency oj p (In the cases of in- 
terest here, ui p will be sufficiently high that this is a rea- 
sonable assumption unless very special precautions are 
taken as in refs. Ull [6^). With Rhf <C Rq, gaussian 
phase fluctuations at oj p are strongly damped; if also 
Rh / Zl , the classical boundary condition at the wire 
ends is approximately just a short, and can be applied us- 
ing the method of images [13] . The resulting current dis- 
tribution for a type II phase slip when I < Xe is shown in 
Fig.EJa) (image phase slips are shown with dashed lines), 
which can be viewed as the adiabatic evolution along 
one period of the lowest energy band C/o($) in fig. ^c) 
(dashed arrow) 28] ; in J J language, the "phase particle" 



(whose "position" is $) jumps from one "potential well" 
of the lowest band Uo{&) to the next [13]. This process 
can only occur, however, if phase fluctuations of the en- 
vironment can excite the system over the barrier. Thus, 
we see that QPS is masked by a low R env (uj p ) Rq 
[13, where for low currents <7 ( < i > ) ~ < i )2 /2ife (dotted 
black line in fig. Etc)) and the wire simply looks like an 
ordinary superconductor with kinetic inductance L k - 

The effect of an external bias current h can be de- 
scribed, just as for JJs, by an additional potential energy: 

U n (I b ,$) = U n ($)-I b $ (25) 

which lowers the potential barrier for phase slip s in one 
direction while raising it in the other 0, [||, HH, HH, H3] 
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[fig. Ete),(f)]. As the barrier is lowered by increasing 
the phase particle has an increasing chance to surmount 
it per unit time due to a (thermal or quantum) phase 
fluctuation. If this occurs, it can either be re-trapped in 
the adjacent potential well by the damping due to Rhf, 
or it can "escape" into the voltage state corresponding 
to a terminal velocity V = 4> (determined by its effective 
mass and the damping) (59j . In the experiments of ref. 
I22I escape rates were measured vs. h and these data 
were fitted to determine an "effective temperature" T e f j 
for the phase fluctuations [ill, H3|. At higher T, it was 
found that T rs T e ff, while at low temperatures T e // 
saturated at a minimum value Tq, and this was taken 
as a signature of QPS [22] (this is very similar to the 
method used to observe macroscopic quantum tunneling 
in JJs [ill). 

In our model, however, we expect an apparent Tq 
which is not directly associated with QPS, but rather 
with the quantum fluctuations of the damped oscillator 
formed by C sh , L k , and R hf [fig. lib)] [gj. Figure Hg) 
shows that this expectation is indeed consistent with the 
observations in ref. [22j for four of the five reported wires. 
We can also compare the average switching current I sw 
observed in ref. [22| with our prediction based on cq. 1251 
Figure. HJh) shows that the agreement is also very good 
for the same four wires. 

Our model also suggests a different explanation for 
another observation in ref. [22j that was was highlighted 
as direct evidence for QPS: the fact that the width of 
the stochastic probability distributions P{I SW ) (obtained 
from many repeated I sw measurements) increased as T 
was lowered. Since the system is overdamped, at high 
T the phase particle moving in the potential Uo(h,$) 
can be thermally excited over a barrier many times (un- 
dergo many phase slips), each time being re-trapped by 
the damping, before it happens to escape into the volt- 
age state. At low T, these excitations are sufficiently rare 
that in a given time it is more likely to experience a single 
fluctuation strong enough to cause escape than it is to ex- 
perience multiple weaker fluctuations which act together 
to cause escape. Just as for J Js, this produces a P{I SW ) 
that broadens as T is lowered [54[ , since fewer phase slips 
are associated with each switching event, and the result- 
ing stochastic fluctuations of I sw are larger. Note that 
in contrast to ref. where these results were explained 
by local heating of the wire by individual quantum phase 
slips, in our model the energy Jf,$o released during a type 
II phase slip is dissipated in the environment Rh / . 

So far we have focused on T <C Tq- However, much 
experimental work has focused on the region near Tq 
where the condensation energy Uc goes to zero. In 
this regime, it turns out that the ratio \e/£, is typi- 
cally less than unity, corresponding to the type I case 
for electric flux penetration (this is in contrast to the 
magnetic case, where the Ginzburg-Landau k — A/£ is 
a temperature- independent material property). The cor- 
responding "type I phase slip" consists of a null in the 
order parameter of dimension £, inside which is contained 
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FIG. 7: Alternative explanation for Giordano's crossover. 
The blue line shows Giordano's fit using the LAMH theory, 
for his wire having A =41nm, and I =80/im from ref. 03 . 
The red line shows our prediction for the same wire, in the 
\e > I, El 3> Es limit. The temperature dependence of all 
quantities was taken from the supplementary information of 
ref. [l^. The vertical dashed lines and filled arrows indicates 
the edges of the two relevant regimes in which the two ex- 
pressions are expected to be valid. This figure compares very 
favorably with fig. 1 of [l3l. In contrast to the conventional 
explanation of the resistive "tail" in terms of a crossover from 
thermal to quantum phase fluctuations, here we explain it in 
terms of a crossover in the temperature dependence of the 
thermal phase-slip barrier only. Note that we have not dis- 
cussed here the detailed behavior in the crossover between 
these two regimes. 



a pulse of electric field with total area (in x, t space) <&o. 
Screening currents j of dimension \e < £, flow around 
this region of nonzero electric field 52j . These events are 
in fact none other than LAMH phase slips 0, H E3] , and 
in most cases a continuous transition occurs from type I 
to type II as T is lowered [61|. 

This picture suggests a new interpretation of Gior- 
dano's original experiments fl3j |. in which LAMH-type 
behavior near Tq was observed to cross over at lower tem- 
peratures to a weaker scaling with T—Tc- At low enough 
T where the type II limit is reached, if Xe > I, the energy 
barrier in our model is w E L /4 (for Es <C E L ) [fig.[5fc)]. 
This scales as oc T — T c (compared to oc (T — T c ) 3 / 2 in 
the LAMH regime), producing the crossover shown in fig. 
[7] (corresponding to the 41-nm wire in fig. 1 of ref. 1 131 ) 
[52j. We therefore interpret the crossover observed by 
Giordano not as a transition from thermal to quantum 
fluctuations, but instead as a crossover only to a different 
temperature scaling of the energy barrier for purely ther- 
mal fluctuations [68j . This picture can also explain the 
mysterious fact that Giordano's anomalous low-T resis- 
tance was not observed in refs. US even though the 
wires had cross-sectional areas nearly 50 times smaller 
and coherence lengths 10 times smaller (which is difficult 
to explain in any theory where the low-T resistance is 
directly identified with QPS). For the wires of ref. [H, 
our calculated stays comparable to £, and never be- 
comes much longer than I, even as T — ¥ 0. Therefore, 



10 




Wire width [nm] 



Wire width [nm] 



FIG. 8: Quantum phase slip energy, (a) Calculated E$x = 
Qo/C±\e, which is the phase slip energy Es for a wire with 
length I — \e (plotted on the right axis for C± =50 pF/m) 
[52| . For a given width, the Es value on the left axis corre- 
sponds to the wire length I on the right axis. In addition to 
9 nm Moo.79Geo.21 H|], 10 nm Al 1j], and 20 nm Ti [TJ, 
several other materials are shown which have low superfluid 
density [J which corresponds directly to low Uc' 12.5 nm 
thick a-Nbo.15Sio.85 [H; 5 nm TiN 0]; 20 nm a-InO* [|; 35- 
nm a-InOx [40(|: 10 nm thick a-Nbo.45Sio.55 HJ. Note that 
our model assumed 2R < £ (not the case for some of the ma- 
terials over the whole plotted range), and treated the wire 
as cylindrical, (b) and (c) show Es vs. wire width for the 
materials and wire lengths of refs. [H and [H[ (sample A), 
respectively, with the vertical dashed line indicating the ex- 
perimental width, and the filled circles the observed values. 
The dashed black line is the Giordano estimate [13l. ll5l . l22i . l44|| . 
the red dot-dash line is the prediction of GZ [3, and the solid 
blue line is our prediction 



no crossover at lower temperatures is expected, precisely 
because the QPS is much stronger in these wires. 

Our discussion has revealed that in a low-Z environ- 
ment, even at T = 0, QPS can be observed in transport 
measurements only indirectly. As long as the avoided 
crossings in fig-HIc) are large enough to suppress Landau- 
Zener transitions, the magnitude of Es has little if an 
effect on the result. However, Es is perhaps the quantity 
of greatest interest, both fundamentally and for possible 
devices. Figure HJa) shows our calculation, for several 
materials, of the quantity Es\ — QqEcXe which is the 
phase slip energy for a wire with I = Xe (we take this as a 
measure of the maximum possible Es in lumped-element 
PS J circuits). We find that Es can in fact be quite large 
for reasonable wire dimensions. This might appear coun- 
terintuitive given the apparent difficulty of observing it, 
but in fact our model suggests that this difficulty has 
been a consequence only of the low-Z environments (at 
uip) used in the relevant experiments. 

Very recently, however, two experimental groups have 
circumvented this problem, and appear to have observed 
Es directly. In the first, using a-InO films, Astafiev and 
co-workers demonstrated the phase-slip qubit of ref. [iil . 
in which the nanowire forms a closed loop. Although 



this can be viewed as Z env = 0, the fixed phase bound- 
ary condition imposed by flux threading a closed super- 
conducting loop allows the PSJ to be biased exactly on 
the avoided crossing of width Es shown in fig. EJc). 
In this case direct spectroscopic measurement yielded 
Es/h ~ 5-10 GHz [40]. The second experiment used 
a-NbSi nanowires embedded in an environment formed 
by highly-resistive Cr nanowires. Although the parasitic 
reactance of these wires likely reduces the impedance at 
high frequencies [54| , the authors were able to observe the 
Coulomb blockade associated with Vc ~ 20 — 100/xV [4lT |. 
Figures [5Jb) and (c) show a comparison of our predic- 
tions with these results, in both cases showing reasonable 
agreement. For comparison we also show the predictions 
of Giordano's model [H, [H [H, H , and that of GZ [yj, 
which are significantly different from both observations. 



IV. DESTRUCTION OF 
SUPERCONDUCTIVITY IN ID 

We now turn to the destruction of superconductivity 
down to T = for short wires with R n > Rq. Previous 
theories have predicted insulating or metallic behavior as 
the wire diameter [14j |. the characteristic impedance Z L 
[l3,[23|, or an external shunt resistor [13] is tuned through 
a critical value, under various conditions. However, none 
can obviously explain a T = transition at R n ~ Rq. In 
all of these theories the predicted transition relies on the 
presence of a form of dissipation which somehow remains 
even as T — > 0, such as anomalous excited quasiparticles 
31 1 , a resistive shunt [27| , continuum plasmon modes 
3113], or the quantum phase-slips themselves [29j . 

We propose an alternative view, in which a T = 
SIT is driven by disorder-induced quantum phase fluc- 
tuations. This is analogous to the SIT observed in some 
quasi-2D systems with low superfluid density p, 4| when 
the sheet resistance Rq > Rq [62]. This 2D disorder- 
induced SIT has been interpreted using the "dirty bo- 
son" model of Fisher and co-workers || , in which disor- 
der nucleates (virtual) unbound vortex-antivortex pairs 
(VAPs), with sufficient strength that these unpaired 
vortices themselves form a Bose-condensate, destroying 
long-range phase coherence and producing a gapped in- 
sulator [5]. This is closely related to the Berezinskii- 
Kosterlitz-Thouless (BKT) vortex-unbinding transition 
in the classical 2D XY model [631 ] . 

To connect these ideas to our system, we recall that 
the imaginary-time path-integral of eq. [TJ] represents an 
isomorphism between closed Feynman paths describing 
the quantum ground state of a ID system at T = and 
statistical configurations of a classical 2D system at fi- 
nite T (5f| [72j • This type of mapping has been exploited 
to predict in ID a T = quantum analog of the classi- 
cal BKT transition [3, 127l l55ll. a phenomenon known as 
instanton condensation [55l |72|] (referred to as "prolifer- 
ation" of phase slips in refs.M E3, 111, III) [Hj]. We can 
apply this idea directly to our type II phase-slips [fill ]: 
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as the disorder is increased, amplitudes for broken phase 
slip-anti phase slip pairs (closed paths in eq. [15]) increase 
in the ground state, until at some critical point they over- 
lap sufficiently to form a condensate. In the dirty boson 
model, the T = critical point at Rq ~ Rq = <I>o/(2e) 2 
corresponds to approximately one vortex crossing for ev- 
ery Cooper pair crossing 0. In our ID case, the corre- 
sponding critical point would appear to be R n ~ Rq . In 
fact, in a recent experiment using "honeycomb" bismuth 
films, consisting essentially of a 2D network of nanowires, 
a SIT was observed precisely when R n of each nanowire 
passed through Rq (by tuning the film thickness) [73$ . 
supporting our proposed connection between the ID and 
2D phenomena. 

This connection may also be able to explain why the 
SIT tended not to be observed in longer MoGe nanowires 
@, [lH . Since the logarithmic interaction required for the 
BKT-like quantum transition only persists to separations 
p ~ Xe [c.f. eq. [M], we might expect to see a disappear- 
ance of the SIT as the wire becomes significantly longer 
than Xe [z3|- Our theory predicts Xe ~ 100-300 nm for 
MoGe wires 5-10 nm wide, which is close the observed 
length range above which the SIT is no longer clearly 
seen. Finally, we note that the observed reduction in Tq 
near the SIT [23|, [2|| , which is not predicted by the dirty 
boson model, may be explained by the gap suppression 
observed in 2D for similar MoGe films [7l|, believed to 
be due to an enhanced Coulomb interaction [75|, [76[ . 



of the phase only. In a low-impedance environment, qua- 
siclassical phase-slip excitations with a normal core arise 
out of this theory (in the same manner that Bloch oscil- 
lations in JJs arise from Josephson tunneling) which are 
electric analogs of the penetration of single magnetic flux 
quanta in 2D. In terms of these excitations, we can pro- 
vide a consistent and in some cases quantitative descrip- 
tion of many phenomena observed for quasi-lD supercon- 
ducting wires in low-Z environments, includin g L AMH 
phase slips [1, Giordano's "resistive tails" |13j. and 
anomalous IV characteristics of ultranarrow MoGe wires 
[TBI I22I I23I |2(| . We also predict values for the phase slip 
energy E$ in good agreement with recent observations. 
Finally, our model suggests a mechanism for the observed 
SIT in short wires when R n > Rq , in terms of a disorder- 
driven quantum phase transition [5], [42| ■ 

An exciting result of our theory is the prediction that 
it should be feasible to achieve large enough E$ to 
build dual circuits to classical Josephson devices, such 
as a quantum standard of current [32l . HH, |36[ dual to 
the Josephson voltage standard, or sensitive electrom- 
eters based on the circuit of ref. |4l] (dual to the DC 
SQUID). Another interesting possibility would be the 
dual of RSFQ digital circuits: a voltage-state logic in 
which Cooper pairs are shuttled between islands (52[. It 
would have no static power dissipation, and could be 
amenable to integration with charge-based memory el- 
ements. 



V. CONCLUSION 

We have described a new theory for quantum phase 
fluctuations in ID built on the hypothesis that flux- 
charge duality 32] holds at the microscopic level between 
these fluctuations and the charge fluctuations associated 
with Josephson tunneling. While many previous theo- 
ries have treated QPS essentially as a quantum version 
of LAMH phase slips (having a normal core inside which 
the order parameter goes to zero [3, [2?], [29|, [3l| ) , in our 
theory QPS is a tunneling process involving fluctuation 
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I. DETAILS ON THE DERIVATION OF AN 
EXPRESSION FOR THE ELECTRIC FIELD 
USING DUALITY 

Here we describe the reasoning behind eqs. 10 and 
11 in the text which allow us to calculate Ec using the 
probability flux density for the macroscopic fluxon wave- 
function. Following ref. 7, we write zero-energy, time- 
independent Schrodinger equations for i\)q and In 
eq. [1] below, Ui = U cp n s can be viewed as the potential 



4> 2 
2£a~ 



Uid 



^ Q (x)=0, 0<x<d 



$ = -i$ dd x 
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v / 2d£ A U I 



$ dim [^(x)d x il) Q {x)\ . 
J = — p-, — t^ttt: = Jc sin < 
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To solve eqs. [T] we use the fluxoid and quasicharge oper- 
ators of eq. [2] For the JJ case, using the semiclassical re- 
sult A = 2m/(2e) 2 n s Q, we obtain: k = y/2(2m)U cp /h, 
which is simply the decay length for a particle of mass 
2m to penetrate into a potential barrier U cp at zero in- 
cident energy. Just as current density J can be written 
in terms of probability flux density for the Cooper pair 
wavefunction, we can also write the fluxon current den- 
sity E in terms of the fluxon wavefunction. The use of the 
ID wavefunction u = i/rip® and the extra factor of 2tt 
for the PSJ in eq. |4] arise from the fact that fiuxons are 
treated mathematically as one-dimensional objects Q, 
whose probability flux is only conserved when integrated 
around the azimuthal angle. 



II. TYPE I PHASE SLIPS 

Figure Q] illustrates a type I phase slip as discussed in 
the text. The type I limit where Xe <C £ can be viewed 
in the following way: is the length of wire whose 
kinetic series capacitance Ck/^E equals its shunt capac- 



energy density for the superconducting state to penetrate 
into the insulator, and Uc (the condensation energy) as 
the potential energy density for the insulating state to 
penetrate into the superconductor. We have also defined 
the specific inductance and capacitance: Ca = A<i and 
Ci = eA. The former can be viewed as the kinetic in- 
ductance x area for the superconducting state inside the 
junction barrier, and the latter the capacitance x length 
of the insulating state inside the superconductor. 
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itance C±Xe- In a type I or type II phase slip, This is 
effectively the length of wire necessary to hold the charge 
that flows onto and off of the core's series capacitance 
Cfc£ during a phase slip (passage of a fluxon through the 
wire). In the type I limit, this charging length becomes 
negligible and drops out of the problem (except insofar 
as the charging of the core boundaries is still required to 
drive the dynamics of the phase evolution). In this limit 
we find that a picture emerges which is different from 
the original intuition of Mooij and co-workers [3| that 
QPS involves tunneling through the same barrier that 
the wire passes over classically during an LAMH (type 
I) phase slip. Instead, in our theory, the type I phase 
slip is real penetration of a single fluxon through an ob- 
servable, localized null in the order parameter, whereas 
a quantum phase slip is coherent (virtual) tunneling of 
a single fluxon through the entire wire, without suppres- 
sion of the order parameter. 



We can use this picture of type I phase slips to derive 
an intuitive result for the so-called "attempt rate" f2 that 
appears in the LAMH phase-slip resistance p], 0, : 
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R ps = Rq 



Ml 



exp 



AF 



(5) 



Recall that the motion of Cooper pairs according to Lon- 
don's hrst equation can be described in terms a classical 
density n s of charges 2e moving with velocity v and mass 
2m, using the relations: A = 2m/n s (2e) 2 and J = n s 2ev. 
Based on the picture above, we treat type I phase slips 
as classical, moving fluxons with kinetic energy described 
by Q 2 /2eA [c.f. eq. 10 in the text], a one-dimensional 
density = l/£, and mass — tc^o/Ci such that 
e = m^/n^&Q and E = ^^n^Vcf,/ R. It is easily veri- 
fied that with these choices, the kinetic energy density 
then reduces to the classical ID electrostatic energy den- 
sity: n ( j > in ( pv 2 1 /2 = (E 2 /2e) ■ A. We can also rewrite 
uifp = {eB 2 ) ■ with = $ /-R£, such that the ef- 
fective mass of a single fluxon is exactly what we would 
classically associate with a flux $o in the (f> direction con- 
tained inside a length £ of the wire. We can estimate the 
thermal activation rate of these fluxons over an energy 
barrier AF for a wire of length I using the simple classical 
expression: 



„ n<hV<hl 

r = ^ r ex P 



-AF 
k R T 



(6) 



To obtain a thermal average for v^, we use: ksT/2 = 
m r j ) (v 2 ^/2, which results in: 



£$0 



(7) 



where = jC^ is the root mean square thermal 

voltage at temperature T across C% = eA/£, the series 
capacitance of a length £ of the wire. Equation [7] ex- 
presses the intuitive conclusion that the average attempt 
rate is simply the incoherent sum of thermal capacitive 
phase fluctuations across Z/£ segments of wire, each of 
length £. Note that this does not include the distributed 
shunt capacitance Cj_ which is associated with fields out- 
side the wire. 




FIG. 1: Type I phase slips. The superconducting gap is sup- 
pressed in an area with dimension £ (length ~ £ in the x 
direction and duration ~ £/c s in time), inside which a region 
of finite electric field exists. Screening "currents" j of width 
~ Ae < 5 now around this region of nonzero flux, just as 
in the 2D magnetic case. This process is none other than an 
LAMH phase slip. Panels (a)-(e) show snapshots in time of 
the evolution. First, a normal region of spatial size £ begins 
to appear (a); then, current begins to flow, charging up Ck on 
the spatial boundaries of this region (b). This produces a lo- 
cal electric field which results in a winding (or unwinding) of 
the phase inside the fluctuation region, until the LAMH sad- 
dle point is reached [J, 0] and the current passes through zero 
(c); The current goes through zero and reverses, discharging 
Ck on the boundaries (d), until the system returns towards 
the uniform state (e). The net result is a single $o "slip" in 
the fluxoid across the wire. Note that unlike type II phase 
slips (Fig. 4 in the text), all of the dynamics occurs within a 
length £. 



a surface energy of the interface between normal regions 
where the flux penetrates and superconducting regions 
where it is excluded (positive for type I, and negative for 
type II). 

In our present 1+1D system, the most direct analog 
to the magnetic case just discussed would be an "electric 
flux" bias. The example of fig. 1(c) in the text, where an 
electric field is produced by magnetic flux in a moving 
frame, illustrates this. The additional term in the Eu- 
clidean action (eq. 15 in the text) associated with work 
done by the source can be written as: 



III. DETAILS ON THE ANALOG TO GIBBS 
FREE ENERGY 

For a superconductor subjected to an external mag- 
netic field H e , the thermodynamics of the field penetra- 
tion is governed by a Gibbs free energy G = F — VH e ■ B, 
where F is the Hclmholtz free energy and B is the actual 
magnetic flux density [7] . The Gibbs free energy accounts 
for the fact that the external flux source must do work 
if the flux is excluded from some or all of the supercon- 
ductor (since it is effectively energizing a smaller induc- 
tance). It is used to describe the different type of flux 
penetration in type I and type II materials in terms of 



S ext = -J (^-j ^da (8) 

and can be understood as arising from the mechanical 
work required to maintain the motion of the magnetic 
source. Analogous arguments can be made about the 
surface energy of normal and superconducting regions: 
for example, in the time dimension, for type I electric flux 
penetration (LAMH phase slips), the order parameter is 
first suppressed before any currents flow, such that the 
source does not provide any energy; for type II phase 
slips, a fluxoid of $0/2 has already appeared across the 
source before the gap begins to decrease near the core. 
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Wire 


A l sw lc DC Lk £jL 
[nm J l/^AJ [KJ [ml J |nrij [1-tlz] [liriz] 


SI 


74 2.37 3.9 67 1.1 2.9 170 


S2 


86 1.4 3.8 66 1.7 1.9 140 


S3 


130 1.42 3.2 61 0.73 4.5 7.0 


S4 


92 0.91 2.9 58 2.2 1.5 180 


S5 


150 4.9 4.6 73 0.51 6.3 0.48 



TABLE I: MoGe wire parameters used in figs. 5(e)-(f). For 
all wires we use p n = 1-8 /j,£l-m from ref. |6|, £ = 4.5 nm from 
ref. [2^ and t = 25eo and a — 0.21 which give a reasonable fit 
to the data in fig. 5(f) of the text. For fig. 5(e) of the text we 
obtain R^f ~ 115Q, from the linear fit shown, with C s h = 1 
fF. 



As the gap goes to zero and the phase evolution along the 
wire is concentrated into a region of length £ [figs. 4(b)- 
(c) ], the effective inductance of the wire (seen by the 
source) is changing while the source is providing nonzero 
such that the source provides energy. Analogs to the 
usual results for flux penetration follow from this, for 
example: a lower critical electric field below which an 
electric Meissner effect would exclude E ea:t , and above 
which an Abrikosov lattice of type II phase slips would 
form [7J. 

Experimentally, the wires are biased with an exter- 
nal voltage or current source, rather than electric "flux" . 
This case can be described by: 



j da 



(9) 



where 3 e xt = + Vs/Z^y with I s and V s the current 
through and voltage across the external source, respec- 
tively. Equation [9] describes work done by the source 
when either: (i) the effective inductance of the wire 
changes while the source is supplying a given current; 
or (ii) when the effective series capacitance of the wire 
changes while the source is supplying a given voltage. 
This expression reduces to that used in ref. [n] when the 
contribution from V s can be neglected, and for a fixed 
source current. 



IV. DETAILS ON CALCULATIONS AND 
PARAMETER VALUES 

Table U lists the parameters used in our calculations 
for wires Sl-5 in ref. (fig. 5(g) and (h) in the text). 
These values are taken from the supplementary material 



of ref. |6|, except for the coherence length £, which we 
take from ref. 0- The values listed for £ in ref. @ were 
obtained only indirectly through fits to LAMH phase slip 
rates, and varied rather widely, from 5-12 nm. 



To obtain the inductive energy El, we calculate the 
Ginzburg-Landau (GL) penetration depth A from the 
sheet resistance using the dirty-limit relation: 



A 



k R T( 



B-LC 



(10) 



where A = fioX 2 , p n is the normal-state resistivity, and a 
is a numerical constant Q . The inductive energy is given 
by: 



E L = 



$2 



HA 



2L k 2/i A 2 / 



(11) 



To calculate the condensation energy Uc, we use the 
value for A and the GL coherence length £ in table U 
along with the GL result for the thermodynamic critical 
field: 



B n = 



27rv / 2A£ 



(12) 



We then have Uc — B c /2^jl . This value is used to cal- 
culate E c : 



En = 



2U C 



e 2tt k^R[K {k^R) 2 + I (k^R) 2 } 



2U C 



exp(— 2«;$i?), 1 



(13) 



(14) 



where Kq and Iq are the modified Bessel functions, and 
is given in eq. [3J Note that although eq. Q3] comes 
from the exact solution to the Schrodinger equation at 
zero energy for all k<s>R, the use of that equation is itself 
only an approximation which may break down for small 
Just as in the macroscopic quantum model of a JJ, 
the use of a Schrodinger-like equation inside the barrier 
is based on the approximation that interactions between 
fluxons (Cooper pairs for a JJ) can be neglected in that 
region. This approximation also underlies the description 
of the phase-slip and Josephson potentials [l(| : 
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Ujj = -Ej cos (f> = 



2 



(\m + l}(m\ + |m - l)(m|) 



(15) 
(16) 



where on the right only terms which transfer one fluxon 
or Cooper pair through the junction are included. 

The calculated quantum temperature shown in fig. 
5(g) of the text is associated with the zero-point fluc- 
tuations of the resonant circuit formed by the kinetic 
inductance of the wire and the lumped shunt capaci- 
tance of the environment C s hi an d is given by [l2j]: 



hjjr 



fkOn 



(2(jJqt) 2 20JqT 
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(17) 



Material 


t £ T c Rn e B c 
[nm] [nm] [K] [fi] [e ] [mT] 


Refs 


Nbo.isSio.ss 


12.5 58 0.21 1400 110 0.48 


[ILiS] 


Nbo.45Sio.55 


10 5 1.2 550 110 24 


[9, 18, 19] 


TiN 


5 10 2 4200 30 7.8 


[20, 21] 


InO* 


20 6 1.5 3500 20 6.2 


[22-24] 


InO* 


35 10 2.7 1600 20 5.4 


[24, 25] 


In 


42 52 4.2 2.9 20 28 


f4] 


Moo.79Geo.21 


8.6 4.5 3.9 210 25 67 


[6, 26] 



TABLE II: Material parameters used in Fig. 



where r = RhfC s h and luq = 1/ ' \JLkC s h- The approx- 
imate equality holds in the overdamped limit u> t <C 1, 
which is well satisfied for the wires of ref. for Rhj ~ 
100ft and C sh -1-10 fF. 

To calculate the predicted switching currents I sw 
shown in fig. 5(h) of the text, we start with the classi- 
cal energy U(h, $ e ) = Uo($ e ) ~ h&o of the lowest band 
shown in fig. 5(c), for a fixed fluxoid $ e across the wire, 
and a fixed current bias lb. The calculation of L/ (<I>e) 
is formally equivalent to the finding the energies of the 
phase-slip qubit @, where in that case the phase bias 
is supplied via an external flux $ e through the closed 
loop. To calculate the energy levels, we diagonalize the 
Hamiltonian: 



H($ e ) = E L 



$0 



-Escos I (18) 



where [Q, $] = ih. We use a discrete fourier grid repre- 
sentation in fluxon states (q\m) = exp(imq) which yields 
[7 n ($ e ) numerically. An example of the resulting poten- 
tial energy for wire SI as a function of $ e is shown in 
figs. 5(e)and (f) of the text. Quantum fluctuations of $ 
at zero temperature about the classical minimum-energy 
equilibrium value, given by eq. [17] above, are also shown. 
We take the average switching current at which the phase 
particle escapes the potential well to be that at which the 
potential well depth is equal to Tq. This condition is used 
to obtain the predicted I sw vales shown in fig. 5(h) of 
the text. 

Table [TT] shows the parameter values used for the plots 
in fig. 6 of the text. For the Giordano 0] and GZ 
results in Figs. 6(b) and (a), we used: 




exp 



-0.3^) il-M 



E 



, ;2 = l8 l_RQkBTc 



Rq 
Rt 



(20) 



where in both cases we have taken unknown numerical 
factors of order unity to be 1. 



V. CONNECTION TO JJ ARRAYS 

The model we have described for QPS, in which the 
quantum fluctuations that allow fluxon tunneling arise 
from an effective mass for the phase which is associated 
with the permittivity e, is conceptually linked to earlier 
work on ID JJ arrays as well as the work of Glazman and 
co-workers in which QPS in a wire was modeled using a 
ID J J array (27j . in fact, one might view our model as a 
continuum version of a ID JJ array, in the sense that such 
an array consists of alternating domains of purely real 
and purely imaginary permittivity (the insulating barri- 
ers, and superconducting islands, respectively) whereas 
we combine the two into a continuous medium with com- 
plex permittivity. This connection can be illustrated by 
re-casting our result for the phase-slip energy Eg in the 
following form (using eq. [T2J : 



E s 



l -yj2V L V c exp 



V L 



2L f 



V c 



2V L 
V c 

Ql 

2C 5 



(21) 
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where = A£ /A and C% = eA/t; are the kinetic in- 
ductance and series capacitance, respectively, of a fluc- 
tuation volume t t = £ ■ A. The quantity Vl is the ID 
analog of the so-called "supcrfluid stiffness" , defined in a 
bulk superconductor for = £ • £ • £ [16[ , which describes 
the potential energy cost (the spring constant) for phase 
fluctuations. We then associate Vc with a kinetic en- 
ergy (effective mass) of those fluctuations. Equation |2"T1 
is nearly identical to the result for the width of the lowest 
quasicharge band in a ID J J array (e.g., ref.[13), with Vl 
playing the role of the Josephson energy Ej, Vc playing 
the role of the charging energy Ec, and //£ playing the 
role of the number of junctions in the array. 

The electric penetration depth Xe also has an ana- 
log in the physics of JJ arrays. In that case, there is 
an electric screening length given by: A c = -^/C/Cn (in 
units of the lattice constant of the array) [28[, where C 
is the junction barrier capacitance, and Co is the shunt 
capacitance to ground (or to a low-Z bias source) of each 
superconducting island. This is essentially identical to 
our result for a distributed PS J (eq. 18 in the text), with 
the barrier capacitance suitably replaced with the kinetic 
capacitance. 



VI. ALTERNATIVE EXPLANATION OF 
GIORDANO'S RESULTS 

Figure 7 in the text shows how the present theory 
may explain the resistance "tails" observed by Giordano, 
which have been interpreted as direct evidence for QPS 
The blue line shows our result for the LAMH resis- 
tance, where we have used the same result for the barrier 
height as Giordano (the LAMH expression scaled by a 
factor of 1/4 H) and our result from eq. [7] for the pref- 
actor. For the red line in the plot we have used eq. EJ but 
with an energy barrier of AF = -El/4, and an attempt 
rate of: 



where V = y / /cbT / 'C ' s h is the rms thermal voltage across 
C s h, for which we have used C 8 h ~ ' £ e// ~5 fF (where 
I = 80/im is the wire length and e e ff ~ 5eo is the av- 
erage permittivity of the surrounding media). We have 
used the results of ref. Q for the temperature depen- 
dence of relevant parameters. The dotted lines indicate 
where each result is no longer valid. This plot compares 
favorably to figure 1 in Ref. 



VII. PSJ CIRCUITS 

Figure [5] shows specific examples of flux-charge duality 
applied to more complicated JJ-based circuits. Panels 
(a) and (b) show the duality between a charge qubit and 



the phase-slip qubit of ref. [||. PSJ-based superconduct- 
ing qubits may be of particular interest since flux and 
charge noise will have their roles interchanged relative 
to JJ-based qubits. Since the excited-state lifetimes of 
present-day JJ-based qubits are thought to be limited 
by high-frequency charge noise, exchanging this for high- 
frequency flux noise (which is thought to be much weaker 
29]) should result in much longer lifetimes. Panels (a) 
and (b) also illustrate how polarization charge on the 
nanowire (produced by a nearby gate electrode) is dual 
to magnetic flux through the junction barrier of the JJ. 
Just as a Fraunhofer interference pattern will be observed 
in the magnitude of Ej vs. flux through the junction (due 
to the Aharanov-Bohm effect) Q , the same pattern will 
be observed in the magnitude of Eg vs. charge on the 
nanowire (due to the Aharonov-Casher effect [8j). This 
may be important for the phase-slip qubit since it im- 
plies charge noise on the nanowire would show up as Vc 
noise in the qubit (dual to Ic noise commonly observed 
in JJ-based qubits [3(|). Panels (c)-(f) show two tunable 
superconducting qubits and their dual circuits. Just as a 
DC SQUID can be used to implement a flux-tunable com- 
posite JJ, the series combination of two PSJs as shown 
can be used to implement a charge-tunable composite 
PSJ. Note that (d) is essentially a tunable version of the 
phase-slip oscillator of Ref. [3l|, and (f) is a tunable 
version of the phase-slip qubit [9j. 

In addition to qubits, where well-defined, long-lived 
energy eigenstates are required in which quantum zero- 
point fluctuations must be kept undisturbed by the en- 
vironment, the circuits shown in (g)-(l) are intended to 
function in a regime where either fluxoid (for JJs) or 
quasicharge (for PSJs) is a classical variable (i.e. where 
quantum fluctuations are small). A well-defined fluxoid 
requires a low environmental impedance at the Joseph- 
son plasma frequency, which is readily obtained using 
resistively shunted Josephson junctions. A well-defined 
quasicharge requires a high environmental impedance 
(^> Rq) at the phase-slip plasma frequency, which is 
much more difficult to realize. In ref. 19, highly-resistive 
Cr nanowires were used to bias the device; in ref. |32| . 
frustrated DC SQUID arrays in an insulating state were 
used. Panel (h) shows the "quantum phase slip transis- 
tor" QPST, first suggested in ref. [H, and implemented 
in ref. 191 This device is an electrometer, dual to the 
DC SQUID amplifier shown in (g). The QPST is simi- 
lar to a single Cooper-pair transistor (SCPT) [34|; how- 
ever, it could have a much higher sensitivity than an 
SCPT, which is limited by the charging energy of the 
JJs (by how small one can make the junction capaci- 
tance). The QPST is instead limited by the kinetic ca- 
pacitance Cfe, whose ultimate limit is the series capaci- 
tance of the wires, which can be much smaller. Panel (i) 
is the Josephson voltage standard, and (j) the quantum 
current standard proposed in ref. [lfj . Under microwave 
irradiation, dual features to Shapiro steps would allow 
locking of the incident frequency / to the applied cur- 
rent I according to / = N /2e, where N is the number 
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FIG. 2: Dual circuits for some well-known JJ devices. 



of parallel PSJs. Such a device would have enormous 
impact in electrical metrology, allowing for the first time 
interconnected fundamental standards of voltage, resis- 
tance, and current (35j |. Finally, panel (k) is a Josephson 
transmission line, a basic building block of rapid single 
flux quantum (RSFQ) digital logic; (1) shows the dual to 
this, in which shunt JJs are replaced by series PSJs, flux 
stored in loops is replaced by charge stored on islands, 
and current bias is replaced by voltage bias. Such circuits 
could be of great practical interest, both because unlike 



RSFQ they have no static power dissipation, and also be- 
cause voltage-state logic could be significantly easier to 
integrate with memory elements than flux-state logic. 
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